Golden rule decay versus Lyapunov decay of the quantum Loschmidt echo 
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The overlap of two wave packets evolving in time with slightly different Hamiltonians decays 
exponentially oc e '', for perturbation strengths U greater than the level spacing A. We present 
numerical evidence for a dynamical system that the decay rate 7 is given by the smallest of the 
Lyapunov exponent A of the classical chaotic dynamics and the level broadening U 2 /A that follows 
from the golden rule of quantum mechanics. This implies the range of validity U > VAA for the 
perturbation-strength independent decay rate discovered by Jalabert and Pastawski [Phys. Rev. 
Lett. 86, 2490 (2001)]. 
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The search for classical Lyapunov exponents in quan- 
tum mechanics is a celebrated problem in quantum chaos 
[1|. Motivated by NMR experiments on spin echoes 
[2 1, Jalabert and Pastawski 0] have given analytical ev- 
idence, supported by computer simulations ||, that the 
Lyapunov exponent governs the time dependence of the 
fidelity 



M(t) = \(ip\exp{iHt)exp{-iHot)\ip)\ 



(1) 



with which a wave packet ip can be reconstructed by 
inverting the dynamics with a perturbed Hamiltonian 
H = Hq + Hi. They have called this the problem of 
the "quantum Loschmidt echo". The fidelity M(t) can 
equivalently be interpreted as the decaying overlap of two 
wave functions that start out identically and evolve un- 
der the action of two slightly different Hamiltonians, a 
problem first studied in perturbation theory by Peres || . 
Perturbation theory breaks down once a typical ma- 
trix clement U of Hi connecting different eigenstates of 
Hq becomes greater than the level spacing A. Then 
the eigenstates of H, decomposed into the eigenstates 
of Hq, contain a large number of non- negligible compo- 
nents. The distribution p(E) (local spectral density) of 
these components over energy has a Lorentzian form 



p(E) = 



2tt(£ 2 + T 2 /4) ' 



(2) 



with a spreading width T ~ U 2 /A given by the golden 
rule Up] . A simple calculation in a random-matrix model 
gives an average decay M oc exp(— Ti) governed by the 
same golden rule width. This should be contrasted with 
the exponential decay M oc exp(— At) obtained by Jal- 
abert and Pastawski ||, which is governed by the Lya- 
punov exponent A of the classical chaotic dynamics. 

Since the random-matrix model has by construction an 
infinite Lyapunov exponent, one way to unify both results 
would be to have an exponential decay with a rate set by 
the smallest of T and A. We will in what follows present 
numerical evidence for this scenario, using a dynamical 
system in which we can vary the relative magnitude of T 



and A. There exists a third energy scale, the inverse of 
the Ehrenfest time te, that is smaller than the Lyapunov 
exponent by a factor logarithmic in the system's effective 
Planck constant. In our numerics we do not have enough 
orders of magnitude between 1/te and A to distinguish 
between the two, so that our findings remain somewhat 
inconclusive in this respect. 

Because T cannot become bigger than the band width 
B of Hq (we are interested in the regime Hi < Hq), a 
consequence of a decay M oc exp[— imin (A, T)] is that 
the regime of Lyapunov decay can only be reached with 
increasing U if A is considerably less than B. That would 
exclude typical fully chaotic systems, in which A and B 
are comparable, and set limits of observability of the Lya- 
punov decay. 

The crossover from the golden rule regime to a regime 
with a perturbation-strength independent decay, ob- 
tained here for the Loschmidt echo, should be distin- 
guished from the corresponding crossover in the local 
spectral density p(E), obtained by Cohen and Heller M. 
The Fourier transform of M (t) would be equal to p(E) if 
ip would be an eigenstate of Hq rather than a wave packet. 
The choice of a wave packet instead of an eigenstate does 
not matter in the golden rule regime, but is essential for 
a decay rate given by the Lyapunov exponent. 

The dynamical model that we have studied is the 
kicked top H, with Hamiltonian 



H = (ir/2T)S y + (K/2S)S 



2 E^ 



tit). 



(3) 



It describes a vector spin (magnitude S) that undergoes a 
free precession around the j/-axis perturbed periodically 
(period r) by sudden rotations around the z-axis over an 
angle proportional to S z , The time evolution of a state 
after n periods is given by the n-th power of the Floquet 
operator 



Fq = exp[-i(K/2S)S 2 z ] exp[-i(w/2)S v }. 



(4) 



Depending on the kicking strength K, the classical dy- 
namics is regular, partially chaotic, or fully chaotic. The 



dependence of the Lyapunov exponent A on K is plotted 
in the inset to Fig. 1 (cf. Ref. EQ)). The error bars reflect 
the spread in A in different regions of phase space, in par- 
ticular the presence of islands of stability. For K > 9 the 
error bars vanish because the system has become fully 
chaotic. For the reversed time evolution we introduce 
as a perturbation a periodic rotation of constant angle 
around the cc-axis, slightly delayed with respect to the 
kicks H , 



Hi = <f>S x J2 5(t 



TIT 



e). 



(5) 



The corresponding Floquet operator is F = 
exp(— icj)S x )Fo. We have set h = 1 and in what follows 
we will also set r = 1 for ease of notation. 

Both H and Hq conserve the spin magnitude. We 
choose the initial wave packets as coherent states of the 
spin SU(2) group pl[ , i.e. states which minimize the 
Heisenberg uncertainty in phase space (in our case on 
a sphere of fixed radius) at the effective Planck con- 
stant h e g ~ if? -1 . The corresponding Ehrenfest time 
is r E = X~ 1 hxS [O. We take S — 500 and average 
M(t = n) = \(ip\(W) n F^\ip}\ 2 over 100 initial coherent 
states ip. 
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FIG. 1. Decay of the average fidelity M for the quantum 
kicked top with K — 13.1 and S = 500, as a function of 
the squared rescaled time {(f>i) 2 - The perturbation strengths 
range between (j> — 10~ 7 and 10 -6 . The straight line cor- 
responds to the Gaussian decay (h) valid in the perturbative 
regime. Inset : Numerically computed Lyapunov exponent for 
the classical kicked top as a function of the kicking strength 
K. Dots correspond to averages taken over 10 4 initial condi- 
tions (see Ref. fllOJ]). The error bars reflect different results 
obtained with different initial conditions. The vanishing of 
error bars indicates the disappearance of islands of regular 
dynamics. 

We first show results in the fully chaotic regime K > 9, 
where we choose the initial states randomly over the en- 
tire phase space. The local spectral density p(a) of the 



eigenstates of F (in the basis of the eigenstates of Fo with 
eigenphases a) is plotted for three different </>'s in the in- 
set to Fig. 2. The curves can be fitted by Lorentzians 
from which we extract the spreading width T. (It is given 
up to numerical coefficients by T ~ U 2 /A, U ~ fiVS, 
A ~ 1/i5). The golden rule regime T > A is entered at 
cj) c w 1.7 • 10~ 4 . For </> <C <fi c we are in the perturbative 
regime, where eigenstates of F do not appreciably differ 
from those of Fq and eigenphase differences can be calcu- 
lated in first order perturbation theory. We then expect 
the Gaussian decay 



M ex exp(-£/V) ^lnMoc (<j)ty 



(G) 



This decay is evident in Fig. 1, which shows M as a func- 
tion of (<pt) 2 on a semilogarithmic scale for <f> < 10~ 6 . 
The decay (||) stops when M approaches M^ — 1/2S, 
being the inverse of the dimension of the Hilbert space. 
This saturation reflects the finiteness of the system and 
eventually prevails at long times independently of the 
strength of the perturbation. 
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FIG. 2. Decay of M in the golden rule regime for kicking 
strengths K — 13.1, 17.1, 21.1 as a function of the rescaled 
time (\?t. Perturbation strengths range from 4> — 10 -4 to 
10 -3 . Inset: Local spectral density of states for K — 13.1 
and perturbation strengths = 2.5 • 10" 4 , 5 • 10" 4 , 10" 3 . The 
solid curves are Lorentzian fits, from which the decay rate 
r ~ 0.8i4> 2 S 2 is extracted. The solid line in the main plot 
gives the decay M oc exp(— Tt) with this value of V. 

For <j) > <j) c one enters the golden rule regime, where 
the Lorentzian spreading of eigenstates of F over those 
of Fq results in the exponential decay 



M ex exp(-C/ 2 i/A) ^lnMa <j) 2 t. 



(7) 



The data presented in Fig. 2 clearly confirm the validity 
of the scaling (Q). There is no dependence of M on K 
in this regime of moderate (but non-perturbative) values 
of 4>, i.e. no dependence on the Lyapunov exponent (A 



varies by a factor of 1.4 for the different values of K in 
Fig. 2). 

We cannot satisfy A < V in the fully chaotic regime, 
for the reason mentioned in the introduction: The band 
width B (which is an upper limit for r) is B = 7r/2 (in 
units of 1/t), while A ^ 1 for fully developed chaos in 
the kicked top (see the inset to Fig. I). For this reason, 
when the perturbation strength <f> is further increased, 
the golden rule decay rate saturates at the band width 
— before reaching the Lyapunov exponent. This is shown 
in Fig. 3. There is no trace of a Lyapunov decay in this 
fully chaotic regime. 
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FIG. 3. Decay of M in the golden rule regime without 
rescaling of time, for K = 13.1, 4> — j-10 -3 , (j=l, 1.5, 2, ... 5) 
(solid curves) and K — 21.1, (j> = 3 • 10~ 3 (circles). Dashed 
and dotted lines show exponential decays with Lyapunov ex- 
ponents A = 1.65 and 2.12, corresponding to K — 13.1 and 
21.1, respectively. The decay slope saturates at (j) fa 2.5TCF , 
when F reaches the bandwidth. 



We therefore reduce K to values in the range 2.7 < 
K < 4.2, which allows us to vary the Lyapunov exponent 
over a wider range between 0.22 and 0.72. In this range 
the classical phase space is mixed and we have coexisting 
regular and chaotic trajectories. We choose the initial co- 
herent states in the chaotic region (identified numerically 
through the participation ratio). Because the chaotic re- 
gion still occupies more than 80% of the phase space for 
the smallest value of K considered, nonuniversal effects 
(e.g. nonzero overlap of our initial wavepackets with reg- 
ular eigenfunctions of i*o or F) should be negligible. We 
expect a crossover from the golden rule decay (0) to the 
Lyapunov decay M 



M ~ exp(-At) => In M ex At, 



(8) 



In conclusion, we have presented numerical evidence 
for the existence of three distinct regimes of exponential 
decay of the Loschmidt echo: the perturbative regime 
(3), the golden rule regime (|7j), and the Lyapunov regime 
(3). The perturbation strength independent decay in the 
Lyapunov regime is reached in our simulation if A < T, 
which prevents its occurence for fully developed chaos in 
the model considered here. Our numerics are limited by 
a relatively small window between A and 1/te (a factor 
hiS « 6). It remains to be seen if the Lyapunov de- 
cay can be observed under conditions of fully developed 
chaos and T < A by increasing S so that 1/te becomes 
larger than T. It is noteworthy that for a Lyapunov de- 
cay M ex exp(— At), the saturated fidelity M^ — 1/2S is 
reached at the Ehrcnfest time te (as can also be seen in 
Fig. 4), so that a Lyapunov decay for t < te rules out 
golden rule decay for later times. Similar investigations 
in strongly chaotic systems with small Lyapunov expo- 
nents (like the Bunimovich stadium with short straight 
segments) are highly desirable. 



10" 



M 



10 



10 



10 




FIG. 4. Decay of M in the Lyapunov regime, for 
<j> = 2.1 • 1CT 3 , K = 2.7, 3.3, 3.6, 3.9, 4.2. The time is rescaled 
with the Lyapunov exponent A, ranging from 0.22 to 0.72. 
The straight solid line indicates the decay M oc exp(— At). 
Inset: M for K = 4.2 and different <f> = j • 10" 4 , j =1, 2, 
3, 4, 5, 9, 17, 25. The decay slope saturates at the value 
4> ~ 1.7 • 10~ for which T ~ A, even though V keeps on 
increasing. This demonstrates the decay law M oc exp(— jt 
with 7 = min(r, A). 
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once r exceeds A. This expectation is borne out by our 
numerical simulations, see Fig. 4. 
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